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Abstract 



We search for CP non-conservation in the decays of r leptons produced via 
e + e~ annihilation at \fs ~ 10.6 GeV. The method uses correlated decays of 
pairs of r leptons, each decaying to the ktt v t final state. The search is done 
within the framework of a model with a scalar boson exchange. In an analysis 
of a data sample corresponding to 12.2 million produced r pairs collected with 
the CLEO detector, we find no evidence of violation of CP symmetry. We 
obtain a limit on the imaginary part of the coupling constant parameterizing 
the relative contribution of diagrams that would lead to CP violation to be 
-0.046 < 9(A) < 0.022 at 90% C.L. This result provides a restriction on CP 
non-conservation in the tau lepton decays. As a cross check, we study the 
decay angular distribution and perform a model-independent search for a CP 
violation effect of a scalar exchange in single r — > 7nr°iJ T decays. The limit 
on the imaginary part of the r scalar coupling is —0.033 < 9(A) < 0.089 at 
90% C.L. 
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I. INTRODUCTION 



Violation of the combined symmetry of charge conjugation and parity (CP) has been of long- 
standing interest as a possible source of the matter-antimatter asymmetry [|J in the Universe. 
Efforts to search for CP-violating effects have concentrated so far on the hadronic sector: CP 
violation in strange meson decay has been the subject of intensive investigation since its first ob- 
servation in 1964 (2j. Recent studies of hadronic || as well as semileptonic kaon decays provide 
precision measurements of the CP violation parameters. Searches for corresponding asymmetries in 
B meson decays are the focus of several large ongoing experiments (7|||. Recent indications of pos- 
sible neutrino oscillations || make it important to re-examine the question of CP non-conservation 
in the leptonic sector. Such violation is forbidden in the Standard Model but appears as a conse- 
quence of various extensions J10| |. Models predicting lepton flavor violation often also predict CP 
violation in lepton decays [JO. 12]. Among theoretically best-known are the multi-Higgs-doublet 
models (MHDM) [13-p"5|]. In this paper, we study CP violation in r decay, in the context of a 
model with scalar boson exchange [J16.17]. The results of the search are general but are easiest to 
interpret for a specific choice of MHDM. Precision studies of muon decay parameters [15,19] show 
no indication for a CP violation in such decay. Previous attempts to study this question in tau 
decays 



2C .21] provide only weak restrictions on the CP violation parameters. 
The search is carried out using data collected with the CLEO detector operating at the Cornell 
Electron Storage Ring (CESR), where r leptons are produced in pairs e + e~ — ► t + t~ . Specifically, 
we study events in which both r's decay to the inr v T final state. In such events, interference 
between the Standard Model process involving W^-boson exchange and a non-Standard Model one 
involving scalar boson exchange can give rise to observable CP-violating terms. For each event 
we compute a quantity for which a non-zero expectation value would constitute evidence for CP 
violation. 

This article is organized as follows. A general description of CP violation in the leptonic sector 
is given in Section O. An observable used to search for CP violation is described in Section III. 
The data analysis and its results are in Sections IV and A model-independent search for scalar- 
mediated r decays is described in Section VI. Derivations of elements of the probability density 
distribution for a pair of r leptons each decaying into Tnr v T final states and theoretical calculations 
used in optimizing the choice of the CP-sensitive variable are given in appendices. 



II. CP VIOLATION IN LEPTON DECAY 



CP violation generates a difference between the partial decay width for a process {i — > f) and 
the corresponding width for its CP-conjugated process [i — > f )■ Any kinematical observable 
£ associated with a decay can be described as a sum of CP-even and CP-odd components: 

£ — £,even £,odd- (1) 

The average value of this observable is equal to an integral of £ over all available phase space 
multiplied by the probability density, P, for the decay. This probability density also can be, in 
general, decomposed into CP-even and CP-odd components: 

P — Peven ~\~ Podd- (2) 
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FIG. 1. Interference between two amplitudes with CP-even and CP-odd relative phases 
5 and 6. 



The average value of £ is: 

< £ >i-*f= J (Seven + Codd)PdLipS = J (£ even + £ odd ) {P e ven + P dd)dLipS 

= J ievenPevendLipS + J CoddPodddLipS . (3) 

Under CP conjugation P odd changes sign and the average value of the observable £ for the decay 
of the CP-conjugated state is: 

< £ >iCP__ f CP= J (Seven + Sodd)P CP ' dhips = J (Seven + Sodd)(Peven ~ Podd)dLipS 

= J £,evenPevendLipS - J £ odd P odd dLipS . (4) 

If P odd is not equal to zero, then < £ >i^f^< £ > i cp^jcp and CP is violated. 



A lepton decay process described by two amplitudes illustrated in Fig. || may have CP-even 
and CP-odd phases 5 and <f) relative to each other. The probability density for such a process is 
given by: 

\A\ 2 = (A x + A 2 e^e i5 )(A 1 + A^e^ 5 ) 

= A\ + A 2 + 2A\A 2 cos (f) cos 5 - 2A\A 2 sin sin 5. (5) 

The last, underlined, term is CP-odd since the phase changes sign under CP conjugation. In this 
example the CP-odd term is not equal to zero if A\, A 2 , sin</>, and sin<5 are not equal to zero. A± 
and A 2 denote the amplitudes and for physical processes they must be different from zero. Thus 
CP-odd term is not equal to zero if the factors e 1 ^ and e lS differ from zero and are complex. We 
discuss in the following section a theoretical model that satisfies these requirements. 



5 



( a ) V ( k ) 



A w = yh^^T^m A H = 

/* A 



r( q, co) 



3350101-002 

*r(P-|) (b) v(k) jt(Pi) 



^°(P 2 ) r(q, £») 




*r u (p 2 ) 



FIG. 2. Amplitude for (a) standard W exchange, (b) scalar exchange. 



A. CP violation in r decays 

A possible scenario for CP violation in r lepton decays is described pl| by the interference of the 
Standard Model decay amplitude mediated by the W boson (amplitude Aw) with the amplitude 
mediated by the charged Higgs boson in the multi-Higgs-doublet modelQ (amplitude Ah)- These 
amplitudes play the roles of A\ and A2 of the previous section. In this scenario, the charged Higgs 
couples to quarks and leptons with complex coupling constants and, thus, there can be a weak 
complex (CP violating) phase (s'mcfi 7^ 0). The usual choice for the CP-even phase 5 is a strong 
phase which arises due to the QCD final state interactions between quarks. In the following, we 
consider only r decays into hadronic final states and a neutrino. The r decay process is described 
by a sum of two amplitudes: Standard Model W exchange illustrated in Fig. §(a) and a scalar 
exchange illustrated in Fig. ||(b). 

To maximize our sensitivity to possible CP-violating effects, we optimize our experimental 
procedures in the context of a specific model that allows us to calculate the matrix element and 
the probability density function for r decay [24|. We have elected to study r decay into the 
tt^tt^Vt- final state due to its large branching fraction, its distinctive experimental signature and 
its relatively simple hadronic dynamics. However, the procedure is general and can be applied to 
a number of other final states. 

The Standard Model amplitude for r decay into two pions via the W exchange can be written 

as: 



A w ~ u(i/)7 M (l-7 5 )u(r) /„ Q", 

\fv\e iS v 



(6) 



where the hadronic current is parameterized by the relative momentum between the charged and 
neutral pions = Pi — p% multiplied by the vector form factor, f v , described by p Breit-Wigner 



x The three-Higgs-doublet model (3HDM) is the least complicated extension of the Standard Model 
allowing for CP violation in r decays |2S 
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shape: 

fv = ZHt , (7 ) 

" s — m 2 + imT(s) ' 

where s is a squared invariant mass of two pions, m n is a pion mass, and m and T(s) are the mass 
and the momentum-dependent width of the resonance, respectively. The latter is defined as: 

r(s) = { %^rf 12 K'XW (g) 
I elsewhere. 



Here, we neglect the contribution from the p' resonance [25]. 

The amplitude for the r decay into vr^vr via a charged Higgs boson in the 3HDM model can 
be written as |il| : 

A H ~ u(y){l - 7s)u(r) -^^[m u Z*X - m d Z*Y] f s M, (9) 

m Higgs 

* . 'IM^ 

A 

where m T , mH iggs , m u , and are the masses of the tau lepton, charged Higgs, u quark and d 
quark, respectively. 

X, Y, and Z are the ratios of the complex Higgs couplings to u, d quarks and leptons relative 
to the Standard Model weak couplings. The overall Higgs coupling to the 7T7r° system is denoted 
by A: 

TYl 

A = 2 ~ [m u Z* X - m d Z*Y]. (10) 

m Higgs 

Because A is complex, the CP-odd phase, (j>, comes from the imaginary part of the coupling 
constant: 

A = K(A) + »9f(A) = | A| (cos + i sin 0). (11) 

We parameterize the hadronic current as a product of a dimensional quantity, M = 1 GeV/c 2 
providing overall normalization, and a scalar form factor f s = \f s \e lSa . Here 5 S is a complex strong 
phase associated with scalar exchange. 

The choice of f s is not unambiguous. We study three cases: one with f s = 1, the second 
with f s described by the ao(980) Breit-Wigner shape and the third with f s described by ao(1450) 
Breit-Wigner shape (see Eq. [?]) with a width given by: 

^ f rn( s - 4 ™\ y/ 2 if S > (2m,,) 2 , 

| elsewhere. 

The matrix element for the r — ► 7T7r°z/ T decay is given by the sum of the W and Higgs exchange 
processes: 

M(t~ -» tt-tt°v t ) = A w + A H , (13) 
and is given by a sum of terms defined by equations |6| and ^: 

M( T - -» tt"tt°u t ) ~ «(A:)7 M (1 - j 5 )u(q, s)\f v \e l8 "Q» + A u(k)(l + 75 )«(g, »)|/.|e M 'M, (14) 
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where q is the 4- momentum vector of the r, s is the polarization of the r and k is the 4-momentum 
vector of the neutrino. 

After calculations detailed in Appendix A the squared matrix element takes the form: 

\M\ 2 T ±^ 0v ~G + s^ (15) 

where the spin-averaged component of the total width is: 

G = 2\f v \ 2 [2{qQ){kQ) - (kq)Q 2 ] + 2|A| 2 |/ a | 2 M0 

+ 4Ste(A)|/ u ||/ a | cos J M T (Qk) - 43(A+)|/„ | \f s \ sin 5 M T (Qk) , (16) 

and the product of the tau polarization s M and the polarimeter vector uj^ describing the spin- 
dependent component of the decay width is: 

= {T2\fv\ 2 M T (2Q^kQ) - k^Q 2 ) ± 2\A\ 2 \ f s \ 2 M T k^ 

T4Me(A)\f v \\f s \ cos5 (Q^kq) - k^qQ)) 
± 49(A+)|/„||/ a | sin^ (Q^kgj-k^qQ)) 

+ 43?e(A)|/„||/ s | sin J e^q a Q^ + 43(A+)|/„| |/ B | cos 5 e fl0l ^q a Q f3 k' y }. (17) 

Here the difference between the strong phases for the vector and scalar exchanges is denoted as 
5 = 5 V — S s . The parameter A + is defined to be equal to A for and equal to its complex 
conjugate, A*, for r + . Underlined terms are CP-odd. 

B. CP violation in correlated decays of r pairs 

For a single r decay we do not know the polarization vector and consequently can not construct 
the spin-dependent term. However, the situation is different for decays of pairs of r's produced in 
e + e~ annihilations, where the parent virtual photon introduces correlations of the t + and r _ spins. 
In the following we denote the momenta of the particles deriving from r + with an additional bar 
symbol in order to distinguish them from the momenta of the r~ decay products. The probability 
density for the reaction e + e~ — * t + t~ can be written as: 

Pt+t-^-^v tt+^p = Gx [(pq) 2 + {pqf + m T (pp)] xG + oj^oj^ (18) 

where p, p, q, and q are momenta of electron, positron, r~, and r + , respectively. The polarimeter 
vectors of the r~ and r + are contracted through the spin-correlation matrix C a ^ (see Appendix B 
for detailed calculations). 

The probability density for r-pair production and the subsequent decay of each r into the 
7T7r° final state consists of all CP-even and odd terms from the single r decay contracted either 
through the spin-averaged production of the r pair (pq) 2 + (pq) 2 + m T (pp) or through the spin- 
correlation matrix C^ u with CP-even and odd terms of the other r decay (Eq. |i~5| ). Due to the 
spin correlation, these contracted terms are no longer proportional to the spins of the individual r 
leptons and, therefore, they do not vanish after integration over all Lorentz-invariant phase space. 
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To design an optimal method for searching for CP violation, we separate the CP-even and 
CP-odd terms in the r production and decay probability density (Eq. |l^). We refer to the sum of 
all CP-even and CP-odd terms as P even and P dd, respectively: 

P = Peven ~\~ Podd- (19) 

Formulae for P even and P dd are given in Appendix C. Both the CP-even and CP-odd parts of the 
cross section are functions of the kinematical quantities that characterize the r + and r~ decays. 

III. EXPERIMENTAL METHOD 
A. Optimal variable 

A common method used in searches for CP violation is to define a CP-odd observable £, such 
as, e.g., a triple product of independent vectors, and then average its distribution over the data 
set. An average of £ different from zero indicates the presence of a CP-odd component, P oc ^(A), of 
the probability density which, in general, depends on scalar coupling A. The choice of the CP-odd 
observable is not unique. One can always add any CP-odd term to any selected CP-odd observable 
to obtain another CP-odd observable. Different observables have, in general, different sensitivity to 
CP violation. To maximize the sensitivity of our search we construct an optimal variable, £, with 



the smallest associated statistical error. Such a variable was proposed by Atwood and Soni [26] 
and by Gunion and Grzadkowski [Effl for other searches. The derivation of £ is described in detail 
in Appendix D. The variable is equal to the ratio of the CP-odd and CP-even parts of the total 
cross section assuming that the absolute value of the coupling A is unity: 



£ = ^ill. (20) 

Peven 

The average value of £ is: 

< £ >= f ^l {Peven + P odd {K))dLips = [ Podd{1 J Podd{A) dLips. (21) 

J "even J "even 

Since P dd is proportional to the imaginary part of A, we can express it as: 

P odd {K) = 3f(A)P odd (l), (22) 

and 

< £ >= Q(A) / Podd{l)2 dLips. (23) 



The integral in Eq. |23j is always larger than or equal to zero, and equality occurs only if the odd part 
of the cross section vanishes. Therefore, we expect that the average value of £ will be proportional 
to the imaginary part of the Higgs coupling constant and will be positive if 9(A) > and negative 
if 9(A) < 0. Monte Carlo simulation of the £ distributions for the three choices of scalar form 
factors and no CP violation^] are shown in Fig. |3](a). The same distributions for the CP violating 
case 9(A) = 1 are shown in Fig. |3|(b). The structure in these distributions is due to the resonant 
structure in the vector and scalar form factors. 



2 We use a modified version of the TAUOLA r decay simulation package [28] to generate Monte 
Carlo samples corresponding to different values of the complex coupling A. 
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FIG. 3. Optimal observable, £, for (a) Monte Carlo with no CP violation and (b) Monte 
Carlo with a maximal CP violation 3(A) = 1. 



B. Experimental issues 



To calculate a value of £ as described in Section III A , we need to know the directions of the 
outgoing neutrinos, or equivalently the directions of the r's. In cases where both r's decay to semi- 
hadronic final states with two missing neutrinos we can use energy and momentum constraints to 



determine these directions with a two-fold ambiguity [29]. We assume that in the e + e~ annihilation 
the two taus are produced back to back and we ignore the effects of initial state radiation. We 
cannot distinguish the correct solution from the incorrect one. We check that the wrong solution 
does not introduce any asymmetry (discussed below) and we sum the £ distributions obtained from 
the two solutions. Finally, we calculate the mean value of the summed £ distribution to search for 
any evidence of asymmetry. 

To assess the possible bias in computed £ values, we use the KORALB/TAUOLA Monte Carlo 
to calculate the mean value of the CP-odd observable for the correct and incorrect solutions sepa- 
rately, as well as for the sum of the two. Within errors, we observe no artificial asymmetry in the 
£ distribution for the false solution in the Standard Model Monte Carlo sample. However, in the 
sample with CP violation, the asymmetry obtained using just the false solution is significantly dif- 
ferent from that obtained using the correct solution. The Monte Carlo procedure used to calibrate 
this effect is described in Section |V|. 



IV. DATA ANALYSIS 

A. Data and Monte Carlo Samples 

The data used in this analysis were collected at the Cornell Electron Storage Ring (CESR) at 
or near the energy of the T(45). The data correspond to a total integrated luminosity of 13.3 fb _1 
and contain 12.2 million t + t~ pairs. Versions of the CLEO detector employed here are described 
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in Refs. |30(] and [|31j]. From this data sample, we select events consistent with e + e~ — ► t + t~ 
interactions where each r decays into the -k-k°v t final state. We refer to such events as p-p events 
since this final state is dominated by production and decay of the p(770) resonance. The event 
selection criteria are summarized in Section IV B| . 

To estimate backgrounds we analyze large samples of Monte Carlo events following the same 
procedures that are applied to the actual CLEO data. The physics of T-pair production and 
decay is modeled by the KORALB/TAUOL A event generator |28[ ], while the detector response is 
handled with a GEANT-based |32j simulation of the CLEO detectors. For backgrounds coming 
from r decays other than r — ► 7nr v T , we analyze a Monte Carlo sample containing 37.2 million 
t + t~ events in which all combinations of r + and r~ decay modes are present, except for our signal 
p-p process. 

Non-r background processes include annihilation into multihadronic final states, namely 
e + e~ — > qq (q = u, d, s, c quarks) and e + e~ — ► T(45) — > BB, as well as production of hadronic 
final states due to two-photon interactions. Backgrounds from multihadronic physics are estimated 
using Monte Carlo samples which are slightly larger than the CLEO data and contain 42.6 million 
qq and 17.3 million BB events, respectively. The background due to two-photon processes is esti- 
mated from Monte Carlo simulation of 37,556 27 —* t + t~ events, using the formalism of Budnev 
et al. All non-tau backgrounds are found to be negligible (see Section IV. C). 

To study the potential bias of the mean value of the observable £ (Eq. |20| ) that can be introduced 
by data selection, we use Monte Carlo samples generated with and without CP violation. The 
"Standard Model" sample (no CP violation) consists of 2.4 million p-p events. The equivalent 
integrated luminosity of this Monte Carlo samples corresponds to four times that of the CLEO 
data set. Each sample of the CP-violating Monte Carlo contains 1.5 million events. 



B. Event Selection. 



Tau leptons are produced in pairs in e + e~ collisions. The event selection follows mostly the 
procedure developed originally for the study of the r — ► pv T decays |54j. At CESR beam energies, 
the decay products of the r + and t~ are well separated in the CLEO detector. We require an 
event to contain exactly two reconstructed charged tracks, separated in ang le by at least 90°. The 
net charge of the event must be equal to zero. Both tracks must be consistent with originating 
from e + e~ interaction region to reject events arising from beam-gas interactions or from cosmic 
muons passing through the detector. Each track must have a momentum smaller than 0.85 Eb eam 
to minimize background from Bhabha scattering and muon pair production. The momenta of all 
charged tracks are corrected for dE/dx energy loss in the beam pipe and tracking system. No 
attempt is made to distinguish from charged leptons or mesons: all charged particles are 
assumed to be pions. 

Clusters of energy deposition in the calorimeter are considered as candidates for photons from 
7T° decay if they are observed in the central part of the detector (| cos 6\ < 0.707, where 6 is the angle 
between the positron beam axis and the photon direction), are not matched to a charged track, and 
have energy greater than 30 MeV. A cluster must have a transverse energy profile consistent with 
expectations for a photon, and it must be at least 30 cm away from the nearest track projection. 
We select events with only four photon candidates. Pairs of photons with a reconstructed invariant 
mass M 77 between —4 and +3 standard deviations (cr 77 ) of the 7r° mass are considered to be 7r° 
candidates. The ir° invariant mass resolution, cr 77 , varies from 4 to 7 MeV/c 2 , depending on the 
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FIG. 4. Energy of the neutral pions: (a) signal Monte Carlo, (b) tau background. 



7T° energy and decay angle. To reject events with spurious 7r° candidates resulting from random 
combinations of low energy clusters, the tt° energy is required to be greater than 250 MeV (see 
Fig. H). Each ir° candidate is associated with the nearest in angle charged track to form a irir 
candidate. Monte Carlo studies show that this method of assignment is correct 99% of the time. 

The event selection criteria successfully rejects background from high multiplicity multihadronic 
events (e + e~ — ► qq). Background from low-multiplicity two-photon (e + e~ — ► e + e~'yy) events are 
rejected by requiring that the missing transverse momentum be larger than 200 MeV/c and that 
the missing mass of the event to be smaller than 7.5 GeV/c 2 (see Fig. ||). Overall, data selection 
efficiency is estimated to be 7.4%. Extensive systematic studies of the background dependence on 
selection criteria cuts have been done in Ref. [Q and confirmed by the present work. 



12 



4 - 



o 
> 

<D 

o 

¥ 2 

> 
LU 



O 

E 

3 
+-> 

C 

<D 

E 
o 



0) 

(0 

1- 
> 

c 

re 



T 1 1 T 



t — i — i — r 



t — i — i — r 



3350101-007 
~\ 1 1 1 



Signal MC 



CLEO II Data 



2.5 5.0 7.5 

MM of the Event (GeV / c 2 ) 



10.0 



FIG. 5. Transverse momentum and missing mass of the event for the data and signal 
Monte Carlo. The vertical line indicates the position of the cut. 
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C. Estimate of the remaining background 



After applying the same selection criteria to the Monte Carlo simulation of e + e~ — > qq and 
e + e~ — ► e + e~77 processes, we estimated that the remaining background from those sources rep- 
resents less than 0.1% of the selected data sample. Consequently we ignore this background in 
our extraction of < £ >. We also find no evidence for contribution of cosmic ray or beam-gas 
interactions in the selected data. The main remaining background is due to the r-pair events in 
which one of the r's decays into pv T while the other decays into t:2tt v t and the photons from one 
of the 7r°'s are not detected. The contamination from this background source is estimated to be 
5.2%. The second largest background contribution of 2.1% is due to one of the r's decaying into 
the K*v T final state producing a neutral pion plus a charged kaon which is mistaken for a pion. All 
other tau decays provide much smaller contributions with the largest being less than 0.7%. The 
total background contamination from tau decays is estimated to be 9.9%. 

V. RESULTS 

A. Calibration of the parameterization after applying selection criteria 

To relate the observed mean value of the optimal observable, < £ >, to the imaginary part of 
the coupling constant A, the 9(A) dependence of < £ > must be known. To simplify the notation, 
in the following we define A to be the imaginary part of the scalar coupling constant A. 

The CP-even terms for each single tau decay are either independent of A or proportional to A 2 , 
while CP-odd terms are either linear functions of A or proportional to A 3 . Therefore, £, which is 
a product of CP-even and CP-odd terms, must have only odd power terms in A expansion and to 
the first order the mean value is proportional to A with a proportionality coefficient c: 

< i >= c x A. (24) 

To check this assertion and to calculate the regions of A where the linear dependence holds we 
generated 21 event-generator-level Monte Carlo samples (no detector simulation) with different 
values of A varying between -1 and 1, and calculated the mean value < £ > for each one. Each 
sample comprised 10 6 events. For each form of the scalar component, the calculated asymmetry 
distribution was then fit to a fifth order polynomial (see Fig. |6|). Using these distribution we find 
the ranges of A for which < £ > is linear in A. As an additional cross-check of our method, we 
search for (and find no evidence of) terms proportional to even powers of A. 

Measurement errors can introduce a bias in the £ distribution (see Eq. ^) as well as change 
the linearity coefficient c (Eq. ^J) due to event selection, reconstruction and resolution. Therefore, 
this coefficient is calculated after applying selection criteria to events generated with full GEANT- 
based detector simulation and pattern recognition software. For each choice of the scalar form 
factor f s we use five samples of 200,000 signal Monte Carlo events generated with different values 
of A in the region where < £ > is linear in A. For each sample, we calculate the average value of 
the optimal observable and plot it as a function of A. For each form of the scalar component, the 
calculated asymmetry distribution is fit to a straight line to obtain the calibration coefficients. To 
check that the selection criteria do not create an artificial asymmetry we calculate the mean value 
of the optimal observable for Standard Model Monte Carlo samples for each choice of the scalar 
form factor. We list these values along with calculated coefficients c in Table |. 
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FIG. 6. A dependence on average value of optimal observable < £ >, for different 
assumptions on scalar form factor f s . 
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TABLE I. Average values of the optimal observable < £ > for the Standard Model Monte 
Carlo and the proportionality coefficient c for the CP asymmetry fits for different scalar form 
factors. 



Form factor, f s 


< £ >, 1(T 3 


c, 


10~ 3 


1 


0.7 ±0.6 


66.! 


3 ±4.3 


BW(a (980)) 


1.0 ±1.1 


586.^ 


I ± 19.4 


BW(a (U50)) 


0.5 ±0.8 


145. 


8 ±7.3 



We observe that the event selection criteria do not introduce an artificial asymmetry in the £ 
distribution. For all three form factors, the mean value of £ for the Standard Model Monte Carlo 
sample is consistent with zero within its statistical error. This fits show no indication that < £ > 
deviates from linearity within the range determined by event-generator-level Monte Carlo study. 
We then use the coefficients from Table | to calculate A. 

B. Observed mean values 

For each choice of the scalar form factor, we obtain a distribution in £, with two entries per event 
for the two solutions for the tau lepton direction as described in Section III. These distributions 
are shown in Fig. 7, with those from the Standard Model Monte Carlo simulations overlaid. From 
these distributions we compute the mean values < £ > after subtracting the average value for 
Standard Model Monte Carlo, which are reported in the first column of Table II. In each case, we 
use the appropriate empirically-determined coefficient given in Table I to derive a value for the 
imaginary part of the Higgs coupling A, as described in the preceding section. These values along 
with the 90% confidence limits on 9(A) are reported in the second and third columns of Table II. 

TABLE II. Average value of the optimal observable £ after subtracting the average value for 
Standard Model Monte Carlo, calculated value of 9(A) and 90% C.L. limits on 9(A). 



Form factor 


< £ >, lO- 3 


9(A), 10~ 2 


9(A), 90% confidence limits 


fs = l 

f s = BW(a (980)) 
f s = BW(a (U50)) 


-0.8 ±1.4 
-0.6 ±2.4 
0.2 ± 1.7 


-1.2 ±2.1 
-0.1 ±0.4 
0.1 ± 1.2 


-0.046 < A < 0.022 
-0.008 < A < 0.006 
-0.019 < A < 0.021 



C. Systematic errors 

In principle, several possible sources of systematic error can contribute to this analysis. We 
found no sizable effect to alter the limits shown in Table [ij. We describe most significant effects 
below. 
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compared to the Standard Model Monte Carlo prediction (solid line) for (a) f s —l, (b) 
f s = BW(a (980)) and (c) f s = W(a (1450)). 
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p (GeV/c) 

FIG. 8. Momentum distributions for (a) 7r _ and 7r + in data, (b) 7i~/tt + ratio in data, 
(c) charged pions, (d) neutral pions. 

1. Detector asymmetry 

The detector can create an artificial CP asymmetry due to the imperfect tracking, detection 
efficiency, and resolution. To check the symmetry of the CLEO detector tracking, an independent 
analysis was done using leptonic tau decays. Such decays have been studied with ~ 0.1% precision 
p5| and at that level show no deviation from the Standard Model. Therefore, any observed 
asymmetry would indicate detector effects. We find the detector response symmetric with respect 
to the charge within a precision of 0.2% and the detector asymmetry, A, in momentum distribution 
is: 

-0.4% < A < 0.2% at 90% C.L. 

2. Track reconstruction efficiency 

We study the systematic effects due to a possible difference of the track reconstruction efficiency 
for 7T~ and tt + as a function of the pion momentum. To estimate the size of this effect we plot in 
Fig. ra(a) the momentum distribution for charged pions in the reaction The ratio 

of these distributions shown in Fig. ||(b) is consistent with 1. The introduction of a ±1<7 slope to 
the ratio plotted in Fig. ||(b) does not change the value of the coefficient c but changes the value 
of 9(A) by ±0.003. We take this as a measure of a systematic error. 

3. Momentum reconstruction 

Another possible source of systematic error is due to an imperfect Monte Carlo simulation of 
the momentum distributions for the tau decay products. Monte Carlo describes the data very well 
[see Fig. ||(c) and Fig. ||(d)]. We estimate the size of the systematic error by introducing artificial 
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±1ct shifts in the slopes of the charged and neutral pions momentum distributions. These shifts 
do not bias the < £ > measurement, but change the calibration coefficient c. The systematic error 
on c is estimated to be 2.2 x 10 -3 and has a negligible effect (~ 10 -5 ) on the value of 3(A). 



4- Background 

Another source of possible distortion of the result is an asymmetry of the £ distribution induced 
by the remaining background. In order to estimate this effect we denote by £5 and the values 
of £ from the signal and background events, respectively. If we denote the number of signal and 
background events by S and B, then, the observed average value of £ is: 

, E is + E tB _ E is S E£ff B 

< « >- s + b ~ ~s~ x JTb + ~b~ x JTb- (25) 

Since the number of background events in this analysis is much smaller than the number of signal 
events, then this equation can be simplified: 

where < £5 > is a mean value due to signal and < £# > is a mean value due to background. We 
estimate < £b > B / S using Monte Carlo: 

< £b > x^ = (0.04 ± 0.04) x 10~ 3 . 

Thus, the systematic uncertainty arising from this source is negligible on the scale of the sensitivity 
of our measurement. 



D. Summary 

Within our experimental precision we observe no significant asymmetry of the optimal variable 
and, therefore, no CP violation in r decays. Due to the uncertainty in the choice of the scalar form 
factor we select the most conservative 90% confidence limits corresponding to / s = 1: 

-0.046 < A < 0.022, at 90% C.L.. 

These limits include the effects of possible systematic errors. 



VI. SEARCH FOR SCALAR-MEDIATED r DECAYS 

The helicity angle, 6^, is defined as the angle between the direction of the charged pion in the 
7T7r° rest frame and the direction of the tttt system in the r rest frame. In Standard Model, the 
helicity angle is expected to have a distribution corresponding to a vector exchange: 

dN 

^ a + bcos 2 e w7r . (27) 



d cosO 



For scalar-mediated decays, there is an additional term proportional to cos 9 nn that corresponds to 
the S-P wave interference and linearly proportional to the scalar coupling constant A. In general, 
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A is complex and the term linear in cos 9 n7T is proportional to the real and imaginary parts of the 
scalar coupling with coefficients c\ and c 2 , respectively: 

^ a + ci5R(A) cos 9^ + c 2 9(A) cos 9^ + 6cos 2 6 n7T . (28) 



d cosO 



The observation of the terms proportional to cosine of the helicity angle would indicate the scalar 
exchange in the tau decays. In the following, we discuss a model independent method used to 
extract confidence limits on the imaginary part of the scalar coupling A. 

In order to calculate the helicity angle we must know the tau rest frame. Due to the unobserved 
neutrino, the tau rest frame can only be reconstructed with a two-fold ambiguity. We can avoid 
such ambiguity by using the pseudo-helicity angle, 9*. This pseudo-helicity angle is obtained by 
replacing the tau rest frame with the laboratory rest frame where it is defined as an angle between 
the direction of tt^ in the irir rest frame and the direction of the tttt system in the lab frame. 
The difference between the pseudo-helicity distributions for the r + and t~ decays is expected to 



have the same form as given by Eq. 28 but with a different numerical coefficients. 

The term including 9(A) changes sign for tau leptons of opposite charges. Therefore, the 
difference of the pseudo-helicity distributions for positive and negative tau leptons has the term 
linear in cos 9* proportional to the imaginary part of the scalar coupling A only: 

dN(r~) dN(r+) 



d cosO* d cos9* 



2c 2 9(A)cos0*. (29) 



The presence of this term indicates CP violation. 

In this study, we use the same data sample as for the previous analysis with the same selection 
criteria except for the requirement of the successful cone reconstruction. The pseudo-helicity 
distribution for t~ and t + is given in Fig. |9](a). 

The structure in Fig. ^(a) is due to the variation in the efficiency as a function of charged pion 
momentum and 7r° energy. To obtain the product of the imaginary part of the scalar coupling A and 
a linearity coefficient c 2 (see Eq. ^9|), we fit the difference of the two pseudo-helicity distributions 
for negative and positive tau leptons to a first order polynomial. To minimize systematic effects 
due to soft pion reconstruction we perform the fit in the region of —0.7 < cos#* < 0.8, which 
corresponds to pions with momentum higher than 0.3 GeV/c. The obtained value of the slope for 
the data distribution is: 

c 2 9(A) = (4.2 ± 3.6) x 10~ 4 

and illustrated in Fig. |9|(b). It is consistent with zero within statistical error. 

To determine the proportionality coefficient c 2 we follow the procedure described in Section 
V. We calculate the slope of the difference of pseudo-helicity distributions after applying selection 
criteria using five samples of 200,000 signal Monte Carlo events generated with different values of 
A. Then, we fit the calculated slope dependence as a function of A using straight line to obtain c 2 : 

c 2 = (107.5 ± 12.6) x 10~ 4 . 

We use this coefficient to obtain the value of the imaginary part of the scalar coupling 9(A): 

9(A) = 0.028 ± 0.037, or - 0.033 < 9(A) < 0.089 at 90% C.L. (30) 

As expected, the limit on the 9(A) is less strict than the one obtained using the optimal observable. 
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FIG. 9. (a) Pseudo-helicity distribution for r~ and r + in data, (b) difference between 
pseudo-helicity distributions for r~ and r + in data. Solid line is a linear fit and dashed lines 
show extrapolation to the region excluded from the fit. 



VII. CONCLUSIONS 

We have discussed a method of searching for the CP violation in the correlated decays of two 
tau leptons each decaying into the -K^ir®v T final state. The limit on the imaginary part of Higgs 
coupling constant is model dependent and depends on the choice of scalar form factor. The most 
conservative choice of scalar form factor f s = l leads to a conservative limit on a scalar coupling: 

-0.046 < 9(A) < 0.022 at 90% C.L. 

Using pseudo-helicity method we obtain a limit on a imaginary part of the scalar component in 
the tau decays: 

-0.033 < 9(A) < 0.089 at 90% C.L. 

Both limits agree with each other and restrict the size of the contribution of multi-Higgs doublet 
model diagrams to the r lepton decay. 

We can relate the limit on 9(A) to the limit on the product of multi-Higgs model coupling 
constants (see Eq. 10): 



-0.026 (GeV/c 2 ) _1 < — ^ — {m u S(Z*Y) - m d %(Z*X)) < 0.012 (GeV/c 2 



m Higgs 
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APPENDIX A: CALCULATION OF THE SQUARED MATRIX ELEMENT FOR 

THE DECAY r~ ir-<n°v T 

The matrix elements and their conjugates have the following forms: 

Af(r" -> vr-vr ^) ~ jyQM Ap(l - 75 )u(g, s) + A/ S u(fc)(l + 76 )«(g > s), (Al) 

AT (r- -> tTttV) ~ Q"u( gj s) 7m (1 - 75 )u(fc) + A*/|«( g , s)(l - J 5 )u(k), (A2) 

M(r+ - tt+ttV) ~ -fvQ»v(q, s) 7m (1 - 75 )u(A;) + A*/ S v(g, - 75 )„(fc), (A3) 

M*(r+ - 7T + 7r°i/ T ) ~ - W«(fc)7M(l " 7&M?> a) + A/£i;(fc)(l + 7 5>(<7, a). (A4) 



The momentum vectors AT, Q, and form factors /„ and f s are defined in section [II A , We used 



the following rules to form complex conjugates: we conjugate all complex numbers, u u, v *-* v, 
and they exchange places; 1 <-> 1, 7 5 <-> — 7 5, 7/J 7yU , 7/j7 5 <-> 7^75- The minus sign at the vector 



hadronic current of the equations (A3) and ( |A4| ) is due to the anti-symmetric behavior of the p 
under the isospin rotation. To form the squared matrix element |M| 2 = M*M, we use: 

u{y)u{y) = k p j p = v(v)v(u), (A5) 



u{t)u(t) = {q a la + m)(l + A/375), (A6) 

w(t)«(t) = (A« - m)(l + A/375)- (A7) 

where s@ is the spin vector(axial) of the tau. Multiplying matrix elements by their complex 
conjugates gives us the absolute value of the squared matrix element: 

\M\ 2 = M*M ~ | jy | 2 W( A« + m)(l + A/37 5 )7„(l - 7s) A P >(1 - 7s) 

+|A| 2 |/s| 2 ( A« + m)(l + A/375)(l - 7s) Ap(1 + 75) 

+A/s/yQ M ( A« + + A/37s)7„(l " 7 5 ) Ap(l + 7 5 ) 

+ A*r s fvQ fl (q a la + m)(l + A/375)(l - 7 5 ) Ap7pXl - 7s), (A8) 



|M| 2 = M*M ~ l/^lV^l^a " "0(1 + A/375)7/ t (l " 75) Ap7,(1 - 75) 
+\M 2 \fs\\q a 7a ~ m)(l + A/375)(l - 75) Ap(1 + 75) 
-A*f* s f v Q»(q a la - m)(l + A/37 5 )(l - 75) Ap7^(l " 7s) 

- MsfvQ»{q a l a ~ m)(l + ^7/375)7p(1 - 75) A P U + 7s)- (A9) 
Reducing to the traces, we can write the matrix element squared as follows: 

|M 2 | ~ 2\f v \\2{q .Q)(k-Q)-(q- k)Q 2 ) + 2|A| 2 |/s| 2 (<7 ■ k) + 4K(A/ 5 / y )m(Q • k) 
+2\f v \ 2 m(2(s • Q){k • Q) - (a ■ fe)Q 2 ) - 2|A| 2 |/ S | 2 m(s • fc) 
+4»(A&#)((« ■ Q){q-k) - (s ■ k)(q ■ Q)) 
+ 4Q(Af s ft)ep afMp sPq a Q»kP, (AW) 
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|M 2 | ~ 2\f v \ 2 (2(q ■ Q)(k -Q)-(q- k)Q 2 ) + 2|A| 2 |/ 5 | 2 (g • k) 
+m(Af s ft)m(Q ■ k) 

-2|/ y | 2 m(2(s • Q)(k ■ Q) - (s ■ k)Q 2 ) + 2|A| 2 |/ 5 | 2 m(s • k) 
-m(Af s fms ■ Q)(q ■ k) - (s ■ k)(q ■ Q)) 
+ ^(Afsfv)^ P s p q a Q^kP. (All) 

Combining these two equations, dropping the overall factor of two and using A + = A for r~ and 
A + = A* for t + gives us the following expression for the squared matrix element of 

\M 2 \ ~ \f V \ 2 (2(q .Q)(k.Q)-(q- k)Q 2 ) + \K\ 2 \f S \\q • k) 
+25ft(A+/s/y )m(Q ■ k) 

±\f v \ 2 m(2(s • Q)(k -Q)-{s- k)Q 2 ) T |A| 2 |/ 5 | 2 m(s • k) 
±25R(A+ f s fv){{s -Q)(q-k) + (s- k)(q • Q)) 
+ 2^{k + f s fr)t^ pS Pq a Q»kP. (A12) 

We can separate the CP-even and odd part in the squared matrix element by using the following 
equivalents: 

^(A+fsfv) = »(A)|/„||/ a | co S 5t»(A)|/,||/ s | sin<5, (A13) 



Q(A + fsfv) = ^)\M\fs\ sm<5±3(A)|/„||/ s | C os5. (A14) 

Now we can rewrite Eq. A12| as: 

\M\ 2 T ±^ „ ~ \f v \ 2 [2(qQ)(kQ) - (kq)Q 2 ] + |A| 2 |/ S | 2 (^) 

2»(A)|/„||/ a | cos 6 M T (Qk) -2$(A+)\f v \\ f s \ sin 6 M T (Qk) 
+s^{±\f v \ 2 M T (2Q^kQ) - k^Q 2 ) T |A| 2 |/ a | 2 M T ^ 
±2»(A)|/„||/ a | cos 5 (Q^kq) - k^(qQ)) 
T23f(A+)|/ v ||/ s | sinJ (Q^(kq) - k^qQ)) 
+2»(A)|/„||/.| sine* e mPl q a Q^ 

+ 29f(A+)|/„||/ a | cos 5 e^q a QP^}. (A15) 



APPENDIX B: r-PAIR PRODUCTION AND SPIN CORRELATIONS 

The matrix element of e + e~ — > t + t~ can be written as: 

M~ [n(r-)7^v(r + )][w(e + )7 / ,u(e")]. (Bl) 

We denote the four- momenta of as q\ and q2, their spins as s\ and S2, respectively, the mass 
of tau lepton as m and the four- momenta of positron as pi and electron as p2- We then write the 
squared matrix element as: 

\M\ 2 ~ Tr{^^ 2ll/ )Tr((l - 75^2)^2 + m)^(l - ^h){h ~ (B2) 
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Calculating the traces we get the following expression: 

\M\ 2 ~ (pi • qi)(p 2 ■ q 2 ) + (pi ■ Q2)(P2 • qi) + m 2 (pi -p 2 ) 

-(si • s 2 )(pi ■ q2){P2 ■ qi) ~ (si • s 2 )(pi ■ qi)(p2 ■ (J2) + (si ■ s 2 )(pi ■ p 2 )(qi • q 2 ) 

-m 2 (si ■ p 2 )(s 2 ■ pi) - m 2 (si ■ pi){s 2 ■ p 2 ) 

~{Pi -P2)(si ■ q 2 )(s2 ■ qi) ~ {Pi -P2)(si ■ giKgg • g2) 

+ (Pi ■ <?i)(<?2 ■ si){p 2 ■ s 2 ) + (p 1 ■ q 2 )(p 2 ■ s 1 ){q 1 ■ s 2 ) 

+{P2 ■ q\){q2 ■ s\)(pi ■ s 2 ) + (p 2 ■ q 2 )(pi ■ si)(qi ■ s 2 ) 

- (<?i ■ q2)(P2 ■ si)(pi • s 2 ) - (qt ■ q 2 ){pi ■ si){p 2 ■ s 2 ). (B3) 

The gauge condition (qi ■ Si) = leads to a vanishing underlined term. Re- writing the squared 
matrix element in terms of spin-averaged production P and spin-dependent part si^s^C^ we 
obtain the explicit form of the spin-correlation matrix: 

|M| 2 ~p + Slti s 2v C^ = 

(Pi ■ <li){P2 ■ q2) + (pi ■ q2)(P2 ■ qi) + m 2 (pi ■ p 2 ) 

+si fl s 2u {-g' MU [(p 1 ■ q 2 )(p 2 ■ qi) + (pi ■ qi)(p 2 ■ q 2 ) - (pi ■ p 2 )(qi • q 2 )] 

- m 2 p^p'i - m 2 p^p v 2 - (pi ■ p 2 )q 2 1 ql + (pi ■ qi)q 2 P 2 
+ (Pi ■ Q2)p 2 qi + (P2 ■ qi)q 2 Pi + (P2 ■ Q2)Px<K 

- (Qi ■ Q2)p 2 Pi ~ (qi ■ q2)PiP2], (B4) 

or 

C» v = -g^[( Pl ■ q 2 )(p 2 ■ qx) + (pi ■ q x )(p 2 ■ q 2 ) - (pi • p 2 )(q x ■ q 2 )] 
-m 2 p^p\ - m 2 p^p v 2 - (px ■ p 2 )q^q\ + (pi ■ qi)q 2 p 2 

+(pi • q2)p 2 ql + (P2 ■ qi)q 2 Pi + (pa ■ ?2)pi9i 

- (<7i ■ Q2)p 2 Pi ~ (qi ■ q2)PiP2- (B5) 
The explicit form of the spin-correlation matrix is equal to: 

C/w = [g^ a - -^q^qa] C Q/3 [gp v - —qpq v ]. (B6) 



APPENDIX C: CP-EVEN AND ODD PARTS OF TAU PAIR CORRELATED 
DECAY RATE INTO (tt-tt ^) - (tt+tt ^) 



The total probability density for the tau pair decays is given by Eq. 18 



P r +T--*ir-ir°v tt+7t°u = G X a X G + OJaC^Uv, (CI) 



where G, G are the spin-averaged matrix elements squared for t~, t + decay (Eq. |l5|), a = (pq) 2 + 
(pq) 2 + m 2 (pp) is the spin-averaged cross section of e + e~ — > t + t _ , uj^ and O u are the polarimeter 
vectors for t~ and r + (Eq. |l^), and is the spin-correlation matrix (Eq. |B6| ). Both the spin- 
averaged matrix element and polarimeter vector of each tau decay contain CP-odd and CP-even 
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terms. Those terms are contracted with the CP-even tau pair production cross section and spin- 
correlation matrix. Therefore, the total CP-odd part of Eq. ^ is a product of the CP-even terms 
of one tau decay contracted with the CP-odd terms of the other tau decay, i.e., it is a linear 



function of the CP-odd terms. Similarly, the CP-even part of Eq. |18| is equal to the CP-even terms 
of one tau decay contracted with the CP-even terms of the other tau plus the the CP-odd terms 
of one tau decay contracted with the CP-odd ones from the other tau decay. This is a sum of 
the non-CP-violating terms of tau decays and the contraction of the CP-violating terms for both 
decays simultaneously, where each term remains CP-even. The matrix element for the tau decay 



is given by Eqs. 15, 16, and 17, where the terms underlined are CP-odd. Thus, the CP-even term 



of the total cross section is equal to: 
p — 

1 even — 

{\fv\ 2 [2(qQ)(kQ) - (kq)Q 2 ] + |A| 2 |/ S | 2 (^) + 2»(A)|/„||/,| cos 8 M T (Qk)} 

XCJX 

{\Iv\ 2 [2(qQ)(kQ) ~ (kq)Q 2 ] + |A| 2 |/ S | 2 (<?ic) + 2»(A)|/ W | |/ s | cos 8 M T {Qk)} 

(-29(A)|/„||/ a | sin <5 M T (Qk)) xax ( -29(A*)|/,,||/ 8 | sin 8 M T (Qk) ) 

{\fv\ 2 M T (2Q^kQ) - k^Q 2 ) - |A| 2 |/ s | 2 M r ^ 
+25R(A)|M|/ S | cos 8 (Q^kq) - k^qQ)) + 2K(A)\f v \\f s \ sin J e fMx p 7 q a Q@k' r } 

■& v - 

{—\fv \ 2 M T (2Q fl (kQ) - k^Q 2 ) + |A| 2 |/ s | 2 M r ^ 
-2»(A)|/„||/ a | cos 8 {Q^kq)-k^qQ)) + 2^(A)\f v \\f s \ sin 8 e vaPl (tQ^} 



+ 



+ 



{ -29(A)|/„||/ a | sin£ (Q^kq)-k a (gQ)) + 2Q(A)\f v \\f s \ cos 8 e^g^F } 

{ 23(A*)|/„||/ 8 | sincT (QjAjcq) - k^qQ)) + 23(A*)|/„||/ fl | cos 8 e^cfQ^ 7 } . (C2) 
Similarly, the expression for CP-odd part of the tau pair production decaying into r Kiftv T state is: 

Podd = 

{\U 2 [2(qQ)(kQ) - {kq)Q 2 } + \K\ 2 \f s \ 2 {qk) + 2»(A)|/„||/.| cos 8 M T (Qk)} 



+ 



+ 



xa x ((-23f(A*)|/„||/ s | sincT M T (Qk))) 



(-29(A)|M|M sin ,5 M T (Qk)) x ax 
{\Iv\ 2 [2{qQ){kQ) - {kq)Q 2 } + \A\ 2 \f s \ 2 (qk) + 2»(A)|/„||/ S | cost! M T (Qk)} 



{\fv\ 2 M T (2Q fl (kQ) - k^Q 2 ) - \A\ 2 \f s \ 2 M T k^ 
+2U(A)\f v \\f s \ cos8 (Q^kq) - k^qQ)) + 2^(A)\f v \\f s \ mi8 e mPl q a Q^} 
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{ 29(A*)|/„||/ 8 | sind (Q^kq)-k^qQ)) + 2S(A*)\f v \\f s \ cos 5 e va ^Q^ } 

+ 

{ -2<3(A)\f v \\f a \ sin S (Q^kq) - k^qQ)) + 23(A*)|/„||/ 8 | cos 5 e^ 7 g Q Q^F } 

{—\fv \ 2 M T (2Q fl (kQ) - k^Q 2 ) + |A| 2 |/ S | 2 M T ^ 
-2K(A)|/„||/ S | cosHQ^kq) -k„(qQ)) +2M(A)\f v \\f s \ sin 5 e ua ^q a Q^}. (C3) 



APPENDIX D: THE OPTIMAL OBSERVABLE 

There is a freedom of choice of the CP sensitive observable £ used in our search.^] To maximize 
the sensitivity of the measurement we need to select £ with a minimal relative statistical error. For 
N events, the statistical error on the average value < £ > is given by: 



A{ _ J <e>-<t» . m 

If the contribution of the CP-odd term is small (i.e., 9(A) is small), the < £ > 2 term is proportional 
to 1 9(A) | 2 (see Eq. E3) and can be neglected. Therefore: 



where P is a probability density of the process and T is a normalization factor, equal to: 

r = J P dLips. (D3) 

The sensitivity S of this method to the presence of a non-zero value of 9(A) is measured by the 
number of standard deviations of < £ > with respect to zero:0 

s= <^> = HPdL EL 

A£ / <g 2 >-<$> 2 ^ ' 

V N 

We can simplify this expression by separating the CP-odd and CP-even parts of the probability 
density P (see Eq. |l9|): 

P = Peven + Podd- (D5) 

The sensitivity becomes: 



3 We assume that the observable is CP-odd. 

4 The average of the observable £ is equal to zero if there is no CP violation. 
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S = VNTR, where R = f g j P ° dd . 

J £ [Peven + Podd)dLips 



(D6) 



We can further simplify R by using the fact that £ 2 P dd is a CP-odd function: 

J i 2 P odd dLips = 0, (D7) 

and 

R = (I jPodd dLips) 2 

J i 2 Peven dLips ' 

The only £ dependent factor in the expression for the sensitivity is R. In order to find an observable 
which is most sensitive to the CP-odd term, we need to find a function £, for which the value of 
R is maximal. We use the functional differentiation method. R is in extremum, when its first 
derivative with respect to £ is equal to zero: 

5R _ 2P odd J £ P odd dLips J £ 2 P even dLips - 2£ P ev en U £. Podd dLips) 2 = Q 



5£ (J e Peven dLips) 2 



or: 



Podd J £ Podd dLips J f P even dLips - £ P even { J £ P od d dLips) 2 = 0. (D10) 

If / £ -Podd dLips = 0, then the condition above is satisfied, but the sensitivity is equal to zero 
(see Eq. |D6| ), thus such case is not a proper solution. If J £ P odd dLips / 0, then we can simplify 
Eq. D10| to obtain: 



Podd J i 1 Peven dLips = £ P even J £ P odd dLips. (Dll) 
One of the possible condition for the extremum of R is: 



iext = A-^, (D12) 
P even 

where A is an arbitrary constant. The constant A must be real for the variable £ to be CP-odd. 
We can decompose any CP-odd observable £ into £, ex t and the reminder £ = £ — £ e xt'- 

£ = + £ ( D13 ) 

±i even. 



Substituting the value of £ into Eq. D£ gives us the following result: 



( f A-^kl- dLips + J £P odd dLipsf 

r = — P 2 — - — : — -7 1 • (Di4) 

J A 2_z£d dLips + 2AJ P odd i dLips + J Peven dLips 

We can simplify this equation using that fact that value of R does not depend on a constant A 



\2 



(see Eq. p!2|) . We define A such that: 

P««| dLipa = 0. (D15) 
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Then, the numerator and one of the terms of the denominator of Eq. D14 vanish. To determine 



the value of A, we multiply left and right sides of Eq. [D12| by P dd an d integrate over the phase 
space: 



/ 



iPodd dLips = A 



p2 

r odd 

P 



dLips + / £P Q dd dLips, 



(D16) 



tPodd dLips 



iPodd dLips - A 



p2 

r odd 

P 



dLips = 0, if : 



(D17) 



A 



J iPodd dLips 



p2 
P odd 

" Peven 



dLips 



(D18) 



Choosing A in this way, we can simplify R, using the property of Eq. D15: 



R 



P 2 

i even 



dLipsY 



p2 
C odd 



dLips + / i 2 P dd dLips 



(D19) 



p2 



PL 



The numerator ( f D ° dd dLips) and the first term of the denominator f T> ° dd dLips are positive 

VJ -i even ' J ^even 

and do not depend on £. R has a maximal value when the last term in the denominator is equal 
to zero. This is possible only when £ is equal to zero. Therefore, the value of R is at maximum 
for £ = A p° dd , and £ is most sensitive to the CP-odd part of the decay rate. Since the sensitivity 
does not depend on a value of the constant A, we choose A = 1 corresponding to: 



Podd 

~p 



(D20) 
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